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In the first part we investigate in detail if the Bose-Einstein condensation scenario for 
magnons can quantitatively explain the observed field- induced magnetic ordering in TlCuCls. 
We use a bond-operator approach to map the spin system onto hard-core bosons and exactly 
account for the hard-core constraint in the dilute limit. We solve the hard-core model within 
the Hartree-Fock-Popov approximation and discuss its validity and the consequences of this 
approximation for the critical properties. In the second part the effects of spin-phonon and 
spin-orbit coupling are discussed within this framework. We show that the experimental mag- 
netization and specific heat data are well described if a certain type of anisotropy is included. 
We also present predictions for the quasiparticle gap which might be tested in the future. 
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1. Introduction 

Antiferromagnetic systems such as integer-spin chains, 
spin-1/2 ladders with an even number of legs or systems 
where the spins form a network of weakly coupled dimers 
share the common property that an excitation gap A be- 
tween the singlet ground state and the triplet state exists. 
An applied magnetic field H leads to a Zeeman splitting 
of the massive triplet with the lowest mode crossing the 
ground state at a critical field Hg = ^jg^B- The ground 
state for H > Hg can be regarded as a Bose-Einstein 
condensate (BEC) of this bosonic mode"'^"'^. This is par- 
ticularly easy to understand for 5 = 1/2 spins where the 
bond-operator (BO) representation 

St2^li±tc.±ti-ie^p^tl^) (1) 

of two neighbouring spins Si^2 m terms of bosonic triplet 
creation (annihilation) operators tj, {ta) is exact Us- 
ing (1) the spin-Hamiltonian can be mapped exactly onto 
a Hamiltonian in terms of the bosonic operators ta ■ The 
magnetic field H becomes the chemical potential fi for 
the magnons and the concept of BEC is in principle di- 
rectly applicable. Note, however, that the magnons are 
subject to the hard-core constraint i^jj^J,, = 0, i.e., only 
one triplet is allowed at the bond i. The magnons there- 
fore form an interacting Bose gas and the problem in 
terms of bosonic operators is in general as complicated 
as that in terms of the original spin operators. The BEC 
concept for the field-induced magnetic ordering in spin- 
gap systems becomes only helpful if the magnon den- 
sity n is small. More precisely, the average distance be- 
tween the magnons I ~ n~^^^ should be much larger 
than the s-wave scattering length a which is the char- 
acteristic length scale representing the influence of the 
repulsive potential. This implies that a/l ~ n^/^a ^ 1 
so that the magnons have to be dilute. In this case the 
well-established gas approximation^ which involves a sys- 
tematic expansion in terms of the small parameter n^^^a 
is applicable and even the finite temperature properties 
of the interacting Bose gas can be studied analytically®. 
In TlCuCls this situation is realized in magnetic fields 



~ 6 — 7 T and this compound has therefore been 
studied extensively in recent years^'^"^^. Inelastic neu- 
tron scattering (INS)^^ has revealed that TlCuCls has an 
excitation gap A w 0.7 meV in zero magnetic field and a 
bandwidth W ^ 6.3 meV. The dimers in this compound 
arc formed by the = 1/2 spins of neighbouring Cu^+ 
ions and weaker interdimer interactions are mediated by 
the Cl~ ions yielding a three dimensional dimer network. 
On the theoretical side it has been shown that the mea- 
sured magnetisation curves can be qualitatively repro- 
duced within the BEC picture'^. In addition, the magnon 
dispersion has been described by using the BO formal- 
ism^. More recently we pointed out that anisotropics in- 
duced by spin-orbit coupling can influence BEC in spin- 
gap systems dramatically. Taking such anisotropics into 
account we have shown that for TlCuCls it is possible to 
obtain good quantitative agreement between the mea- 
sured magnetisation curves and those calculated within 
the BEC framework^^. 

In the present paper we want to extend"'^^ with respect 
to the following aspects: In section 2 we will show how 
the hard-core constraint can be taken into account for 
dilute magnons beyond the mean-field level thus improv- 
ing the results in Ref.^. In particular, this approach al- 
lows us to calculate the magnon-magnon scattering am- 
plitude directly. In section 3 we will discuss the validity 
of the Hartree-Fock-Popov approximation (HFPA) which 
we use to solve the hard-core boson model. We will also 
investigate how this approximation affects the critical 
properties of the model. In section 4 we discuss the in- 
fluence of crystal-field anisotropics for the BEC. We also 
investigate in more detail the case of a staggered g ten- 
sor and/or antisymmetric (Dzyaloshinsky-Moriya) spin 
interactions proposed in Ref.^^ to explain the measured 
magnetisation curves. In particular we will show results 
for different orientations of the magnetic field. In section 
5 we compare the results of our theory for the specific 
heat with experimental data. The last section presents a 
summary and conclusions. 
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2. Bond-operator formalism 

The BO technique starts from the strong coupUng 
ground state \s) where each dimer at bond i forms a 
singlet \i,s). The operators then create local triplet 
excitations \i, a) = \i, s) with \i, +) = — , \i, — ) = 
lii) and \i,0) = (Iti) + UT>)/V2. For TlCuCls an ef- 
fective Heisenberg-typc Hamiltonian has been derived 
in^' containing the intradimcr coupling J and three 
interdimer couplings J^, Ja2c, Jabc- In these works it has 
been argued that the exchange paths corresponding to 
these interdimer couplings arc most important. This has 
been confirmed in a recent electronic structure calcula- 
tion^^. We have checked that including additional ex- 
change paths which might be also of some relevance^^ 
does not significantly improve the results presented here. 
We therefore restrict ourselves to the model considered 
jj^9, 10 Using (1) and retaining only terms bilinear in the 
triplet operators one easily derives^ 

H = ^i^{Ak + agiJ.BH) t^Jka + ^ {tkat-ka + /i.e.) | 

(2) 

- fk + 9k, Bk = fk + gk, fk = Ja cos + 
kz), 9k = 2JabcC0s{ka; + kz/2) cos{ky /2) , 



with Ak = J 
Ja2c cos{2kx 4 



a = —a where —it < kx,ky < tt and — 27r < kz < 27r. 

the 



By the Bogoliubov transform tka = Uktka 
Hamiltonian gets diagonal with eigenvalues 



^ka = \JaI - Bl - a9iiBH 
and Bogoliubov coefBcients 



2 2 
Uk,Vk 



±1 

2 



Ak 
2oJko 



(3) 



(4) 



A good fit of the measured magnon dispersion^ ^ can be 
obtained from (3) using J, Ja, Ja2c, Jabc as fitting param- 
eters. Because these parameters are not very sensitive to 
the gap value whereas the magnetic properties for fields 
H > Hg crucially depend on A we prefer to fix the gap 
and use it as a constraint in the fitting procedure. Re- 
sults for A = 0.65, 0.8 and 1.0 meV are given to the left 
of each column in table I and are consistent with^. 

Next we discuss the renormalization of the dispersion 
(3) due to the hard-core constraint. Following^ this con- 
dition can be taken into account by introducing an infi- 
nite on-site repulsion between the bosons: 



Hu = tlJljjtiatip, U ^ oo . 



(5) 



The corresponding scattering vertex Tap.apiK) where 
K — {k,Lu) is the total momentum and energy of the 
incoming particles can be calculated exactly in the dilute 
limit by a summation of ladder diagrams. This yields 



T{K) 



(6) 



For magnetic fields H > Hg only particles near the band 
minimum at Qq = (0,0, 27r) are excited so that the en- 
ergy and momentum dependent T{K) can be replaced 
by the constant Vq = r(qrQ,0). Because is close to 1, 



vq is approximately given by the magnon bandwidth W. 
However, it is known that the hard-core constraint (5) 
leads to a renormalisation of the Bogoliubov coefficients 
(4) and also to a renormalised triplet spectrum^ 



kO 



S„(fc,0))'-(Bfe + E„ 



(7) 



Here Z^^ = 1 — dY^n/duj is the quasiparticle residue, 
S„(fe,a;) the normal and the anomalous self-energy. 
Values for these quantities and the renormalised superex- 
change parameters are given in brackets in table I. Al- 
though the superexchange parameters are considerably 
renormalised, the shape of the dispersion is only slightly 
changed because Sa and Z are almost momentum 
independent here and, in addition, a fit to the mea- 
sured dispersion is performed. The most important con- 
sequence is a renormalisation of the scattering amplitude 
vq by more than 20%. The dispersions for all 3 gap- values 
are shown in Fig. 1 in comparison to the INS-data. In all 
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Fig. 1. Dispersion relation for TlCuCls measured by INS (cir- 
cles) in comparison to the theoretically calculated ones with 
fixed gap values A = 0.65, 0.8, 1.0 meV. 



cases the experimental data are well described and dif- 
ferences between the 3 fits are only visible close to Qq. 

3. HFP approximation 

As we are interested in temperatures T < A <C Uo we 
can treat vq as temperature independent and it is also 
sufficient to take only the lowest triplet mode [a = +) 
into account. For simplicity we define = Qfeo — The 
Hamiltonian for the lowest triplet mode is then given by 



H 



\ t ''^O 

Mo) ifcifc + y 



^k+q^k' 

i,k',q 



^k+Jk'-a^ktk' (8) 



where = 9i^B{H — Hg). To allow for the symmetry 
breaking in the condensed phase we introduce new op- 
erators tk = Ck + iSk,q„il where 77 is a real number and 
no = 77^ the condensate density. To diagonalise (8) we 
have to deal with the terms which are cubic or quar- 
tic in the new operators Cfc. The simplest way to treat 
these terms is the one-loop approximation which yields 
the Hartree and the Fock diagram for the quartic term 
(see Fig. 2). Therefore this approximation is often called 
Hartree-Fock-Popov approximation (HFFA). As a result 
we find H = Hc + Hun + Hhmn where 



H. = 



VO 2 

Y^o- Mono 



J. Phys. Soc. Jpn. 



Full Paper 



Jesko SiRKER et al. 3 



Table I. Superexchange parameters, self-energies, scattering amplitude fo and effective mass m for dispersion (3) and (7) (in brackets). 



A [meV] 


J [mcV] 


Ja [mcV] 


Ja2c [mcV] 


Jabc. [meV] 


Sn [mcV] 


Sa [mcV] 


Z 


vq [mcV] 


m [1/mcV] 


0.65 


5.52 (1.83) 


-0.24 (-0.30) 


-1.57 (-1.81) 


0.16 (0.58) 


(1.21) 


(0.46) 


1 (0.91) 


8.01 (9.91) 


0.25 (0.26) 


0.8 


5..52 (1.77) 


-0.23 (-0.21)) 


-1.50 (-1.82) 


0.15 (0.57) 


(1.25) 


(0. 15) 


1 (0.91) 


7.77 (9.11) 


0.32 (0.32) 


1.0 


5.51 (4.78) 


-0.21 (-0.27) 


-1.56 (-1.81) 


0.45 (0.55) 


(1.20) 


(0.43) 


1 (0.92) 


7.39 (8.89) 


0.41 (0.42) 




Fig. 2. (a), (b): Hartree and Fock diagram, (c), (d): Additional 
one- loop diagrams due to the terms cubic in c^^\ 

iJlin = i{2nvoiq + vonori- ^lQ7i){c\^- Cq^) (9) 

Huxin = |^fc4cfe - ^ (4clfc + h.c)j I 

with n being the density of non-condensed magnons. 
Here Ak = e/e — /xq -I- Sn with the normal self-energy 

511 = 2vQn -\- 2vono and the anomalous self-energy 

512 = VoUq. By a standard Bogoliiibov transform ifbiiin 
can be diagonalised and = {A^. — £^2)^''^ is the quasi- 
particle spectrum. In addition Hn^ has to vanish. 

It is important to note that the HFPA is parametri- 
cally justified in the present case: The magnon density 
n = n + hq is small and diagrams taken into account 
in the HFPA are the leading diagrams in a systematic 
expansion in n. However, when T — > Tc from below, 
where Tc is the critical temperature of the Bosc gas, it is 
well known that the condensate density no in the HFPA 
does not vanish. Therefore the total density jumps and 
the HFPA falsely predicts a first order phase transition^ . 
This is due to the fact that any perturbative approach 
fails at the critical point where the scale for fluctuations 
diverges. By estimating when the next-leading diagrams 
become as important as the diagrams taken into account 
one finds^ that the HFPA is justified when 

\Tc — T\^ an '/^Te . (10) 

As on^/^ is our small parameter, the temperature re- 
gion around the critical point where this perturbative 
approach is not justified is usually very small. 

In many papers about EEC in TlCuCls much interest 
has focused onto the fact that Hc{T) - Hg ^ T'^' with 
^ ~ 2^ where is the critical field where BEC occurs. 
This has lead to a discussion why the "critical exponent" 
(/) deviates from 3/2 which is the result within the HFPA 
provided that the magnon dispersion is quadratic. In the 
same way one can also consider n^T) T'^ and finds 
from the experimental data^ that this "power law" yields 
a good fit with a similar exponent </> ~ 2. In HFPA one 
finds nc{T) ~ T^/^. Therefore it has been speculated^' 
that the HFPA fails, i.e., that the "critical exponent" is 
renormalised due to quantum fluctuations. 

First, it is worth to note that we are dealing with quan- 
tities which are not dimensionless. Therefore the expo- 
nent is dictated by dimensional analysis. If the disper- 
sion is given by e fc" we will find in rf-dimensions that 
UciT) <^ T^Ajot^ Of course, the proportionality factor is 



not necessarily temperature independent as in the HFPA. 
But it can only depend on the dimensionless quantity 
aT^/"- (or equivalently an^/^ which is small for a dilute 
gas) where the scattering length a, is the only length 
scale in our problem. Corrections of this kind have in- 
deed been found for interacting Bose gases^^. As a con- 
sequence there is no longer a simple power law. However, 
for the BEC of magnons there is also another effect which 
makes a simple power law invalid even within HFPA: The 
dispersion is certainly quadratic for small excitation en- 
ergies but the range of validity of the quadratic approx- 
imation might be rather small. As we pointed out^^'^^ 
this is indeed the case for TlCuCls where the quadratic 
approximation is only justified for T < 1 K which is well 
below the experimental temperature range^. To see this 
we show in Fig. 3 nc(T) calculated within the HFPA 
using the triplet dispersions from Fig. 1. Whereas for 




T[K] 



Fig. 3. riciT) within HFPA with dispersions from Fig. 1. Upper 
inset: For T < 1 K the dispersions arc almost quadratic and 
ric ~ T'^/^ (dashed lines). Lower inset: ndT) is well fitted by 
a "power law" (symbols) with exponents 1.92 (A = 0.65 meV), 
1.88 (A = 0.8 meV) and 1.82 (A = 1.0 meV) for 2 < T < 4 K. 

T < 1 K we indeed see that n,. ^ T^/"^ there is no power 
law at higher temperatures. However, as shown in the 
lower inset of Fig. 3 it is still possible to obtain an ex- 
cellent fit by a "power law" even at higher temperatures 
provided that the considered temperature range is suf- 
ficiently small. For the temperature range in^ we find 
an exponent ~ 1.8—1.9 which is in good agreement 
with experiment and a recent work^* where also a real- 
istic dispersion has been used. We should not take this 
agreement too serious because the HFPA is not justi- 
fied at the critical point. In addition, we will argue in 
the next section that anisotropics play a crucial role and 
that actually no phase transition occurs. What Fig. 3 
nevertheless does show is that in the experimental tem- 
perature range nc{T) does depend on microscopic details 
and no universal power law exists. 
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4. Spin-phonon and spin-orbit coupling 

In-'^^ we have shown that the magnetisation M{T) cal- 
culated with the dispersion (7) and vq ~ 10 meV as ob- 
tained by the BO formalism does not fit the experimen- 
tal data. However, our calculation of vq in section 2 does 
only include magnon-magnon scattering. From Raman 
spectroscopy-'^^ and sound attenuation experiments^^ it 
is known that phonons are important. Because optical 
phonon modes exist at energies comparable with the en- 
ergy scale of the magnetic excitations"'^ it is not clear 
if the effect of magnon-phonon scattering can be calcu- 
lated perturbatively. We have therefore used vq as a fit 
parameter and obtained best agreement with experiment 
for Vq ~ 25 meV as shown in Fig. 4a. A renormalisation 
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Fig. 4. (a) Experimental magnetisation curves (symbols) for H = 
5.9, 6.0, ■ ■ • , 7 T from Ref.'^ compared to the theoretically calcu- 
lated (solid lines) with A = 0.67 meV, dq = 25 meV. (b) The- 
oretical magnetisation curves as in (a) (solid lines) and v^ith an 
additional exchange anisotropy 7 = 0.01 meV (dashed lines). 



of Vq by a factor 2 — 3 might be caused by a reduction 
of the bare magnon bandwidth due to polaronic effects. 
But even with vq — 25 meV we can obtain good agree- 
ment with experiment only for T > Tc- At low-T our 
theory still overestimates M{T) by 50%. ari^^^ ~ 0.1 so 
that HFPA should be justified according to (10) apart 
from a region ^ ±0.5 K around the critical point. The 
approximation therefore cannot explain the failure of our 
theory at low-T. 

In any real magnetic system there is some kind of 
anisotropy which reduces the symmetry. In a system 
without magnetic field and anisotropies we have the 
usual SU{2) symmetry. By a magnetic field this sym- 
metry is reduced to U (1) around the magnetic field axis. 
Spontaneous breaking of t/(l) occurs at the BEG transi- 
tion and is responsible for a gapless Goldstone mode in 
the phase with no ^ 0. Any kind of anisotropy will in 
general break U{1) explicitly so that there is no longer 
a Goldstone mode. However, depending on the type of 
anisotropy there might be still a Z2 symmetry (chang- 
ing the sign of the triplet operator) so that a transition 
between a phase with and without condensed magnons 
is still possible. Finally, if the anisotropy breaks Z2 no 
phase transition will occur and for a small symmetry 
breaking anisotropy we expect the phase transition to be 
smeared out to a crossover. 



Here we want to discuss two kinds of anisotropy. 
First, we want to consider an exchange anisotropy (EA) 
within a dimer, i.e., we consider J ^ { J^, J^, J^}. Using 
triplet operators, performing the Bogoliubov transforma- 
tion with parameters as in Eq. (4) and considering only 
the lowest triplet mode we find 



(11) 



Jy. This kind 



as perturbation to (8) where 7 cx J, 
of perturbation can also originate from a "single-ion 
anisotropy" for the triplets ~ D{Sf + Sl)^ + E[{Sf + 
5*1)^ — {Sf + 5*1)^] where Si^2 denote the spins within 
one dimer. On the other hand consider a Dzyaloshinsky- 
Moriya anisotropy (DMA) ~ D ■ (Si x ^2) within the 
dimer. Transforming this type of interaction into triplet 
operators we find 



2, pert 



(12) 



where the wave vector depends on how the DM- vector 
D varies in space^^. In particular we want to consider the 
case where Qq corresponds to the minimum of the triplet 
dispersion. Both types of anisotropy yield additional con- 
tributions to the Hamiltonian (9) and results for the 
quasiparticle spectra and the additional constraints due 
to the linear term are given in table II. When compared 
to the case without anisotropy we see that the EA yields 
a small shift 7/fo in Uc and also a small quasiparticle 
gap Aqp = \/8yvQriQ in the condensed phase. At the crit- 
ical point, Aqp is zero. Fig. 4b shows that apart from 
the shift in ric the shape of the magnetisation curves is 
basically unaffected. On the other hand a DMA has a 
dramatic effect even if it is tiny because it smears out 
the phase transition to a crossover. As already shown 
in^^ it is possible to obtain excellent agreement with the 
measured magnetisation curves for -ff || 6 when a DMA 
7 = lO"'^ meV is included (see Fig. 5a). In an ideal 
TlGuCla crystal the centre of each dimer is an inver- 
sion centre. 7^0 therefore requires small lattice dis- 
tortions. If these distortions are of such kind that D is 
oriented along a specific direction throughout the crys- 
tal it would be possible to restore U{1) symmetry by 
applying the magnetic field along the same axis. In this 
configuration a sharp phase transition would still occur. 
We therefore compare here also with experimental results 
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Fig. 5. (a): Experimental magnetisation curves as in Fig. 4a and 
theoretically calculated with g = 2.06^^, A = 0.72 meV, vq = 
27 meV and 7 = IQ-^ meV. (b): As in (a) with g = 2.26^^. 
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Table II. Comparison between exchange and DM-type anisotropy 





cxcliaii^'c anisotropy (- > 0) 


DM-type auisotropA- (- ^ 0) 


Hlin = 


-^or; - 277/ + q {2von + Vorio) = 
Two cases: rj = (iion-condciiscd phase) 
and // ^ (coudoiist^d phase) 


-/ior/ - 7 + {2v()n + 'Uo'rao) = 
ahvays eondeiised magnous. r/ 7^ 


Phase transition 




no phase transition 


non-condensed phase 


Ek = y/ (£fe - /xo + 2hvo)'^ - 




condensed phase 


Ek = Viek + 27)2 + 2(efe + 47)t;ono - 





for H ± (102)19 in Fig. 5b. As for H || b the magnetisa- 
tion curves show only a smooth increase at temperatures 
below the minima and no sharp phase transition. There- 
fore a component of D _L ff seems to exist also for this 
configuration. This could be possibly explained by small 
lattice distortions leading to domains with different ori- 
entation of £). In this case a component D J- H could 
exist for each field configuration. The theoretically calcu- 
lated M{T) show reasonable agreement with experiment 
also for H _L (102) if we change the (/-factor from 2.06 
to 2.26 according to ESR^^'^o. This shows that the vari- 
ation in 7 with field direction and also additional crystal 
field anisotropics seems to be relatively minor. Although 
we could certainly improve agreement with experiment 
by inchiding such effects we have not done so because 
the number of fit-parameters would be too large. 

In Fig. 6a we show the condensed and the non- 
condensed density separately for the same parameters as 
in Fig. 5a. The transition from the condensed to the non- 
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Fig. 6. (a) Condensed and non-condensed magnon densities, (b) 
Aqp as function of T for different H and as function of H for 
T = 1 K, 3 K and 5 K (inset). The parameters are as in Fig. 5a. 



condensed phase is smeared to a crossover in a region ^ 
1 — 2 K around the former transition point. Analytically 
we find that the temperature region AT where the den- 
sities with 7 7^ deviate significantly from those without 
such an anisotropy is given by AT/T ~ {'j/vo)^^^/n. 

As the temperature range where quantum fluctuations 
are important is much smaller according to Eq. (10) 
the HFPA is parametrically justified here even in the 
crossover region. This does not mean that corrections to 
the HFPA do not exist but these corrections can be cal- 



culated everywhere in a perturbative manner. 

Finally, we show in Fig. 6b the quasiparticle gap Agp. 
As a function of temperature Agp has a minimum ba- 
sically at the same point where also the magnetisation 
curves in Fig. 5a have their minima. The theoretically 
calculated Agp as a function of magnetic field at fixed 
temperature does qualitatively agree with recent ESR 
measurements^'^. Quantitatively the gap predicted by our 
theory is about a factor 2 smaller than the one measured 
by ESR for || b at H = 9T. Two explanations are pos- 
sible: First, within our theory it is also possible to obtain 
reasonable agreement with the measured magnetisation 
curves for 7 ~ 5 • 10~^ meV by changing vq and the exci- 
tation gap A accordingly. In this case Agp ^ 0.2 meV in 
agreement with ESR. Second and more important, it is 
not clear if the gap measured in ESR is the pure quasi- 
particle gap. Within our theory we would expect the gap 
measured in ESR at fields above the "critical field" (min- 
imum in the gap function) to be given by a combination 
of one- and two-magnon excitations. In this case the ESR 
gap would be a combination of the real quasiparticle gap 
Aqp and 2Agp^^. To avoid these ambiguities INS mea- 
surements of the gap are desirable. 



5. Specific Heat 

To calculate the specific heat Cv{H) we need the en- 
ergy E+{H) of the lowest triplet component {a — +) in 
the presence of a DMA 7 and a magnetic field H. From 
Eq. 9 and table II one easily finds 



E+{H) 



/io"o - {n + no)vo - 2«o"-o - 4:Vonnn (13) 



where 7 only enters by the condition for the vanishing 
of the linear term (see table 11) . As we want to calculate 
the specific heat for temperatures T ~ A we also have 
to take into account the contribution of the triplet mode 
a = 0. At = all modes are of the same importance 
and the DMA yields only tiny corrections which we can 
ignore. Therefore the energy E{0) at zero field of each 
triplet mode a = +, — , is simply given by 

E{0) = -2von^ + ^^iflko+4von)hl (14) 

k 

hi = l/(exp[/3(nfeo + 4t;ori)]-l). 
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As phonons yield a large contribution to the specific heat 
which we do not want to consider here, we only compare 
theoretical and experimental results for Cv{H) — Cy(0). 
In this quantity most phonon contributions should be 
eliminated. As the a = mode is basically unaffected by 
the magnetic field we find 

d[E+{H)-2Em 



Cv{H)-Cv{Q) = 



dT 



(15) 



The results are shown in Fig. 7a in comparison to the 
experimental data^^. The agreement is good. The overes- 
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Fig. 7. (a) Measured specific heat^^ (per copper ion) for H = 
7, 8, • ■ • ,12 T (circles) and theoretically calculated (lines), (b) 
Calculated specific heat for DMA 7 = 5- lO"^,--- ,10-^. In 
both graphs the subsequent curves are shifted by 0.05. 



timation of the peak heights and slight underestimation 
of the widths indicates that 7 seems to be slightly larger 
than assumed here. In addition we also see some devia- 
tions at higher temperatures particularly at higher fields. 
Here we should remember that at such energies we have 
to use in principle the energy and momentum dependent 
scattering amplitude '^{K) instead of the constant vq. 
The dependence of peak height and width on 7 is shown 
in Fig. 7b. By increasing 7 the peak gets smaller and 
broader. Note that without anisotropy there would be a 
singularity in Cv{H) — C'v'(O) at the critical point due 
to the failure of the HFPA. 

6. Conclusions 

To summarise, we have shown how to incorporate 
the hard-core constraint into the BO formalism if the 
magnons are dilute. For TlCuCls we have found that 
the dispersion is only slightly modified when compared 
to^ where the constraint has been completely ignored. 
This is due to the following facts: (a) The quasiparticle 
residue and the self-energies turn out to be almost mo- 
mentum independent and (b) a fit to the measured dis- 
persion is performed. However, even in this case the cor- 
rect treatment of the constraint is important when calcu- 
lating the magnon-magnon scattering amplitude vq. We 
have solved the hard-core boson model using the HFPA 
and have shown that this approximation is valid apart 
from a small region around the critical point. Even if 
TlCuCls is assumed to be a system without anisotropics 
no power law ric ^ can be expected in the exper- 



imental temperature range^ because the quadratic ap- 
proximation for the dispersion works only for T < 1 K. 
Next, we have discussed how tiny EA or DMA influence 
BEG. An EA yields only a small shift in ric and a small 
quasiparticle gap in the condensed phase but leaves the 
shape of the magnetisation curves otherwise unchanged. 
On the other hand a DMA has a dramatic effect and 
smears out the phase transition to a crossover. The dif- 
ferent effects of these anisotropics can be understood by 
the symmetries of the system. With a DMA 7 = 10^^ 
meV we have achieved good agreement with experimen- 
tal data for M{T) and Cv[H, T) if H \\ b. When rescaled 
by the g-factor the agreement with M{T) for H _L (102) 
was also reasonable. This shows that a component of the 
DM-vector D J- H seems to exist for each field config- 
uration and that additional crystal field anisotropics are 
relatively minor. Within our theory we expect a quasi- 
particle gap Agp ^ 0.1 meV for H \\ b, H ^ 10 T which 
might be tested in the future. 
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